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Abstract. We present an analysis of large-scale structure from the Spitzer 
Wide-area Infrared Extragalactic legacy survey, SWIRE. The two-point angu¬ 
lar correlation functions were computed for galaxies detected in the 3.6-micron 
IRAC band, on angular scales up to a degree. Significant evolution in the clus¬ 
tering amplitude was detected, as the median redshift of the samples increases 
from z = 0.2 to 0.6. The galaxy clustering in the GALICS semi-analytic mod¬ 
els was compared with the observed correlation functions and found to disagree 
with the data at faint flux limits. 


1. Introduction 


The large-scale structure of the universe is a result of the gravitational growth of 
dark matter density perturbations. The detailed structure is dependent on both 
the cosmological parameters, which can be constrained by measurements of the 
cosmic microwave background for example, and the evolution of the dark matter 
distribution, as modelled with N-body simulations and semi-analytic techniques. 

It is known that the distribution of galaxies is biased with respect to the dark 
matter, but the details of the bias is poorly constrained. Observations of the 
large-scale structure of galaxies, and its evolution, can be used to understand 
this relationship between galaxy formation and the mass field. _ 

The Sp itzer Wide-area Infrared Extragalactic survey fSWIRE: lLonsdale et al.l 
l2nn.3l. 1200411 is the largest of the Spitzer Space Telescope’s six legacy programmes. 
When completed, we will have observed a total area of 49 square degrees, split 
between six fields, in all seven of Spitzer’s imaging bands (3.6-160 ftm). SWIRE 
is being used to study the history of star formation and the assembly of stellar 
mass in galaxies, the nature and impact of accretion in active galactic nuclei, 
and the influence of environment on these processes. 

The area and depth of SWIRE combine to produce a survey of significant 
cosmological volume. Out to the estimated median redshift of z = 1, the sur¬ 
vey proper volume is 0.04 Gpc^ - larger t han the Sloan D igital Sky Survey 
(DR3 spectroscopic sample, median z = 0.1: lAbazaiian et ahl 2005). SWIRE is 
sensitive to both star-forming and passively-evolving galaxies in the same vol¬ 
ume, extending out to redsh i fts of 2-3 and across comoving spatial scales up 
to 100 Mpc. In lOliver et ahl (|2004l l we presented the first detection of galaxy 


clustering with Spitzer, measuring a two-po int angular correla tion function at 
3.6 /xm from our validation data. Similarly, iFanEf et alJ 1 2004h presented their 
IRAC angular correlation functions from the Spitzer First Look Survey. Here 
we begin to extend our work to the larger SWIRE survey. 
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Figure 1. (a) Two-point angular correlation functions at 3.6 gm for three 

flux limited samples within the SWIRE ELAIS-Nl field with S' 3.6 > 32 /iJy, 
80 /rJy and 200 gJy. The best-fitting power-law models, w{9) = ® — C, 

are shown for each sample, (b) The amplitude, A, of the angular correlation 
function as a functio n of K-h&nd magn itude, compared with the evolving 
models (light lines) of iRoche et al.l ll2003tl . The heavy/dotted line is the e = 0 
model normalized to the SWIRE data (filled circles). 


2. The two-point angular correlation functions 


The restframe emission of galaxies in the H- and it'-bands is dominated by 
older stellar populations and traces the stellar mass of galaxies. At the median 
redshift of the SWIRE survey, this corresponds to the observed 3.6 fim band, 
thus a sample of galaxies selected at 3.6 /xm is essentially selected on the basis 
of the total mass in stars. 

The ELAIS (European Large-Area ISO [Infrared Space Observatory] Sur¬ 
vey) N1 field has an area of ~9 sq. deg. and ha s been observed in all seven 
Spitzer bands. The Version 1.0 source catalogue (|Surace et al JliH) has been 
used to select three flux-limited samples at 3.6 /xm. A faint sample was selected 
such that the flux limit of Ss.g > 32 /xJy was well-above the survey detection 
limit, ensuring a high source reliability (signal-to-noise > 40) and a uniform 
selection function; this sample covers an area of 1 sq. deg. A medium sample 
with S' 3,6 > 80 /xJy covers 2 sq. deg. and a bright sample has S' 3.6 > 200 /xJy over 
8 sq. deg. Note that the faint and medium samp les are conta i ned w ithin the 
bright sample and are not independent. Eollowing lOliver et al.l ( 2004l i. sources 
were removed from the initial catalogue if they were identified as point sources 
in the Two-Micron All Sky Survey (2MASS) or had the 3.6/4.5-/xm colours of 
stars. The remaining stellar contamination was estimated to be < 3%. 

The two-point angular correlation function, w(9 ), was calculated f or each 
of the three samples following the same techniques as lOliver et al. |(| 200 l) , using 
the Landay-Szalay estimator. The correlation functions are plotted in figure la. 
A power-law model, w{6) = A9^~'^ — C, was fitted to the data over 9 > 0.003 deg 
( 11 '') - the 6 " photometry aperture limits the reliability of the data on smaller 
scales. With a fixed slope of 7 = 1.8, the best-fitting correlation amplitudes. A, 
were (1.8±0.2) x lO'^ at 32 /xJy, (5.8±0.4) x lO'^ at 80 /xJy and (13.9±1.0) x lO-^ 
at 200 /xJy. 
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Figure 2. The SWIRE 3.6-/im angular correlation functions compared with 
the GALICS simulations for the three flux-limited samples. The shaded regions 
are the l-cr error bounds on w(&) from the mock catalogues; the data points 
and best-fitting curves correspond to figure la. 


3. Comparison with GALICS mock catalogues 


GALICS (Galaxies In Cosmological Simulations) is a hybrid model of galaxy evo¬ 
lution which combines high-resolution N-body simulations of the dark matter 
content of the universe with semi-analytic p rescriptions to desc ribe the fate of 
the baryons within the dark matter halos (jHatton et al.ll?003l i. The simula¬ 
tions have 256^ particles in a 140 h^Q^ Mpc box, with a minimum halo mass of 
2 X 10^^ Mq. Within each halo, some fraction of the gas mass is cooled and 
turned into stars which then evolve. The spectral energy distributions of these 
model galaxies are computed by summing the contribution of all the stars they 
contain, tracking their age and metallicity. A mock catalogue is generated by 
projecting a cone through the simulation at a series of timesteps (redshifts), and 
calculating the properties of the galaxies ‘observed’ in the cone. The GALICS 
project have made available several of these cones, from which we have con¬ 
structed mock catalogues of the SWIRE survey, in l-sq.-deg. patches. 

We have calculated the two-point angular correlation functions for the GAL¬ 
ICS catalogues, using three flux-limited samples corresponding to the SWIRE 
data described in the previous section. In figure 2 we compare the GALICS cor¬ 
relation functions (shaded regions) with the SWIRE results (points; figure la). 
At the brighter flux limits (Ss.e > 80 /uJy & 200 /uJy) the clustering in the 
simulations is consistent with the data. For a fixed 7 = 1.8, the amplitudes of 
the best-htting power-laws differ by < 2 -cj between the data and the simulations 
in each sample. The larger uncertainties in the model correlation functions are 
due to the small size (1 sq. deg.) of the simulations compared with the data 
(8 sq. deg. at 200 /uJy) - there are significantly fewer sources contributing to the 
estimation of w(d). In the faintest sample, 53.6 > 32 ^Jy, the model predicts 
stronger clustering than we observe on scales > 1 arcmin. The amplitude of the 
GALICS correlation function, A = (2.7 ± 0.2) x 10“^, differs by > 3 -cj from that 
of the SWIRE data. 
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4. Clustering evolution 


The angular correlation function is the projection along the line of sight of the 
spatial correlation function, and is dependent on both the redshift distribu¬ 
tion and luminosity function of galaxies in the survey and on the evolution of 
the spatial clustering. At fainter flux limits the survey probes to higher redshifts 
and larger volumes, and this reduces the strength of the projected clustering. 
This is shown in figure lb where we plot the amplitude, A, of the angular cor- 
relation function a gainst the limiting magnitude for a range of A'-band surveys 
( Roche et al.ll2nr)^ . We used the GALICS simulations to estimate the equivalent 
AT-band limit for our 3.6-/rm-selected samples, and plot these SWIRE clustering 
amplitudes on the same figure. Our new data have much smaller errors in A 
and are in agreement with the AT-band surveys. 

Also plotted in figure lb (light lines) are three models of clustering evolution, 
^(r, z) = (r/ro)“^(l+^)“*'^^^^ with e = 0, —0.4 & —1.2, from iRoche et al.l ( 

We note that these models appear to be offset from the data for K < 18, and we 
have renormalized the e = 0 model to agree with the SWIRE amplitudes (heavy 
line). This renormalized model is consistent with most of the existing K < 18 
data and, significantly, a linear interpolation of the model passes through the 
2MASS data points at AT < 14 (dotted line). 

We make several conclusions based on these results. Eirst, the correlation 
length ro (or the amplitude A) of stellar mass selected samples (at AT-band &: 
3.6 //m) appe ars to be smaller t han that of the /-band selected sample used 
in the models (jRoche et al.ll2nn^ . Second, the scatter in the ground-based AT- 
band data limits their ability to discriminate between the models (figure lb). 
The SWIRE flux limit, = 3.7 //Jy, corresponds to an equivalent K ~ 20.5 
and we will soon measure the correlation function to this depth, with smaller 
errors than the existing data. Third, the observed amplitude of the correlation 
function at faint flux limits is inconsistent with the GALICS catalog ues (figure 2) 
and a ppears to be more complex than the power-law models of iRoche et al.l 
( 2003t l in figure lb. Detailed modelling of the redshift distribution, luminosity 
function and clustering evolution is being developed, which will enable us to 
gain further understanding of the observed galaxy clustering. 
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